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THEORETICAL AND EXPERIMENTAL ANALYSIS OF ONE-DIMENSIONAL COMPRESSIBLE 
FLOW IN A ROTATING RADIAL- INLET IMPELLER CHANNEL 
By Seymour Lieblein 
SUMMARY 

An analysis of the one -dimensional compressible flow in a rotating 
radial-plane impeller channel was conducted in order to provide an 
insight into the characteristics of the passage mean flow under the 
influence of centrifugal forces and losses. From a theoretical inves- 
tigation of the flow in an illustrative impeller channel with convergent- 
divergent area variation, the behavior of the flow along the channel was 
found to be generally similar in trend to the flow in a stationary 
convergent-divergent nozzle. The critical (sonic) section of the rotat- 
ing channel occurred upstream of the geometric throat. The effect of 
the losses on the flow was similar to the effect of a reduction of the 
flow area. 

A further one -dimensional analysis was conducted of the flow 
in an experimental radial-inlet impeller containing static-pressure 
taps along the stationary front shroud. The behavior of the mean 
flow along the impeller passage was generally similar to that of 
the flow along a rotating radial channel in which the effective flow 
area in the inlet region varied with the operating point. At neg- 
ative angles of attack, sonic velocity was attained in the inlet 
region of the impeller where a throat had formed because of the sep- 
aration of the flow. Calculated maxi mum weight flows based on con- 
ditions at the critical radius compared favorably with experimental 
maximum weight flows over the wheel-speed range of the impeller. 
Impeller-inlet losses were found to be large at the higher weight flows. 


INTRODUCTION 

In an attempt to understand the fundamental-flow behavior of cen- 
trifugal compressors, several theoretical and experimental investiga- 
tions of the flow within centrifugal-impeller passages have recently 
been conducted (references 1 to 6) . Considerable attention has been 
given to the two-dimensional impeller types with radial inlet and dis- 
charge, in particular. These studies have included experimental meas- 
urements of over -all performance and passage flow (references 3 to 5) 
and theoretical calculations of passage and blade-surface flow distri- 
butions (references 1, 2, and 6) . 
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One of the characteristics of these impeller analyses is the 
assumption that the flow proceeds without loss and completely fills the 
geometric area of the passage at all times. Experimental results have 
indicated, however, that, particularly in the inlet regions of impellers, 
friction and mixing losses and boundary-layer separation exert signi- 
ficant effects on the behavior of the channel flow. 

As a preliminary approach to the incorporation of such "real-flow" 
effects in an impeller analysis, the usefulness of one -dimensional con- 
siderations in studying the actual passage mean flow of impeller chan- 
nels was investigated at MCA Lewis laboratory. A one -dimensional 
analysis is helpful for evaluating the various types of losses, for 
locating regions of large loss, and as an initial basis for the intro- 
duction of friction and loss terms in theoretical developments. For 
methods of analysis such as in references 4 and 6, the use of more 
realistic passage mean flows could result in more accurate determina- 
tions of velocity distributions along blade surfaces. 

An academic analysis of the ideal flow in a rotating radial-plane 
channel whose area variation was selected to represent the effective 
flow area existing in experimental impellers is presented in the first 
part of this report. Based on an illustrative impeller channel, the 
effects of centrifugal force and losses on the general one-dimensional 
flow variation and on such factors as location of the critical section, 
maximum mass flow, and shock formation axe investigated. 

In the second part of this report, the one-dimensional analysis is 
applied in conjunction with impeller front-shroud static -pressure 
measurements to an investigation of the flow in the passages of a 
radial-inlet impeller whose over -all performance was reported in ref- 
erence 5. The general variations of the flow are discussed, and cal- 
culation of the effective flow area in the impeller -inlet region is 
made. The choking of the impeller flow and variations of impeller 
losses are analyzed. 


SYMBOLS 


The following symbols are used herein: 


A area 


a 


velocity of sound 


specific heat at constant pressure 


E energy number 
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F energy loss or work done in overcoming losses 
g acceleration due to gravity 
J mechanical equivalent of heat 
K constant 

M flow Mach number 

impeller -tip Mach number (tip speed/a Q ) 
n polytropic exponent 

P absolute total pressure 
p static pressure 

Q heat energy added 

R gas constant 

r impeller radius 

T absolute total temperature 

t static temperature 

V relative velocity 

W mass flow 

w weight flow 

a angle of attack 

y ratio of specific heats 

6 ratio of inlet total pressure to standard sea-level pressure 
Tj ad adiabatic efficiency 

Q ratio of inlet total temperature to standard sea- level temperature 
p mass density 

\|r flow function 


k 
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angular velocity 


Subscripts: 


1 impeller inlet 


c 


critical (sonic) condition 


m maximum 


o 


inlet total (stagnation) conditions 


T impeller tip 


t throat 


The prime refers to irreversible (loss) conditions. 


THEORETICAL ANALYSIS 


Basic Equations 


In a one -dimensional theory of the steady flow of a perfect gas in 
a channel, the state of the fluid is considered to be sufficiently 
defined by average values of its components such that each flow variable 
is a function of a single parameter or coordinate of position. The 
conventional variables - pressure p, temperature t, density p, 
and relative velocity V - as functions of a single space parameter in 
steady flow are connected by four relations derived in reference 7, which 
when altered to include the effect of the centrifugal force field of a 
rotating passage are given by: 

Conservation of energy 



(1) 


Conservation of mass 


d(pVA) = 0 


( 2 ) 


Equation of motion 



( 3 ) 
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Equation of state 


8(p/pRt) = 0 


(4) 


The variations of pressure, relative velocity, and relative Mach number 
along the radius of a channel are then obtained from equations ( 1) to 


The variations of velocity, pressure, and Mach number are thus seen to 
be functions of the individual contributions of the gradients of wheel 
potential energy, friction loss, heat transfer, and flow area with the 
square of the Mach number as the principal parameter. 

For isentropic flow along a rotating impeller channel (dF = 0, 

8-Q = 0) , the equations of motion and energy become identical. The 
ideal state of flow at any radius is then obtained by integrating equa- 
tion (3) isentropically between the inlet of the channel and a point in 
the channel to yield 


If there is no prerotation at the inlet to the rotating passage, the 
velocity terms on the right side of the equation are equal to the 
kinetic energy of the absolute velocity immediately before the entrance 
to the rotating channel, so that 


(4) (appendix A) as the logarithmic differential equations 





r C 2 + (r-l)M 2 J dF - (r-l)(l + rM 2 ) dQ> +[2 + (r-l)M 2 ] ^ 


(7) 


V 2 , r P 
2 + r -i p 



r Pp oi 2 r| 


( 8 ) 


r-i p x 2 


v 2 1 r p 
2 r-i p 


= V*i 


2 


(9) 
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Making use of equation (4) and the isentropic relations, 



Equation (9) can he written in terms of stagnation conditions upstream 
of the channel inlet as 


r-i 



where Mj is defined as impeller-tip Mach number and is equal to 
aar T /a 0 . The ideal state of flow can then be calculated at any point 

along the channel over a range of rotational speeds and weight flows by 
the use of equation (10) in conjunction with the continuity condition. 

The applicability of the general validity of the one -dimensional 
approximation (equations (5) to (7), or (10)) depends on the degree of 
variation of flow conditions across a plane normal to the direction of 
flow. The velocity variations in impeller channels (because of blade 
loading, secondary flows, and boundary layers) are recognized as being 
greater than those normally encountered in tubes and nozzles for which 
the one-dimensional assumption has been found adequate. Calculations of 
energy and efficiency based on full channel area and average velocity 
would therefore contain inherent inaccuracy. However, if the number of 
blades is large and if corrections for effective flow area are made, the 
one -dimensional approximation might yield reliable information concern- 
ing the variation of the average flow along the radius of a rotating 
channel. 


Flow in Illustrative Impeller Channel 

The principal effects of energy transfer and area variation on the 
basic behavior and restriction of the flow in rotating channels can 
best be demonstrated by first investigating the flow of an ideal isen- 
tropic fluid and illustrating the analysis for a representative channel. 
The effects of losses and friction can then be exa min ed as causing 
departures from the ideal configurations. 
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Geometry of channel . - The radial variation of geometric area nor- 
mal to the relative flow path for the illustrative impeller channel was 
selected to represent the effective flow-area variations existing in 
the channels of the radial impellers used in references 3 and 5. The 
geometric area of these impellers was designed to increase along the 
flow path from inlet to outlet. In the inlet region of the impeller, 
the true flow area is not clearly defined because of the variation of 
the magnitude of the relative area of the approaching air as the impel- 
ler weight flow is varied. As shown in figure 1, at low weight flows, 
the relative approach area is less and at high weight flow, the relative 
approach area is greater than the blade -inlet area. A continuous vari- 
ation of area must occur across the inlet, however, as the flow adjusts 
to the channel shape. At the high weight flows, for example, the adjust- 
ment involves a separation of flow from the driving surface, the forma- 
tion of a throat, and then mixing and reattachment of the flow. Evi- 
dence of inlet -area reduction in experimental radial impellers at high 
weight flow is found in references 3 to 5. In order, therefore, to 
employ a fixed-area variation representative of the effective flow area 
of experimental impeller channels, a convergent -divergent configuration 
was chosen for the illustrative channel. For convenience, a parabolic 
variation of area was prescribed. The area variation, equivalent cir- 
cular nozzle and general configuration of the fictitious impeller corre- 
sponding to this area variation with negligible axial curvature, is 
shown in figure 2. Although the effective throat in the experimental 
impeller does not occur as far downstream as indicated in figure 2, the 
exaggerated theoretical location was chosen to accent the effects of 
the throat formation. Blade thickness in figure 2(c) represents both 
actual blade and separated areas. The fluid is assumed to follow the 
mean flow path between blades (direction for infinite number of blades) 
and to completely fill the passage at all times. Inlet, throat, and 
discharge sections of the channel are shown by stations 1, t, and T, 
respectively. Stagnation conditions are assumed at station o. 

Description of flow . - The state of flow along the channel of the 
illustrative impeller (fig. 2) was calculated over a range of rotational 
speeds and mass flows by means of equation (10) and the continuity con- 
dition. The variations with radius ratio of pressure ratio, velocity, 
and Mach number for increasing rates of flow from zero to maximum are 
shown in figure 3 for an impeller-tip Mach number of 1.0. The trends of 
variation of pressure, velocity, and Mach number axe thus seen in gen- 
eral to be si m ilar to those of a stationary convergent -divergent nozzle. 
The instantaneous rise in pressure at the blade inlet is due to the fact 
that the angle of the arbitrarily defined blade shape at the inlet does 
not correspond to the angles of the entering air at the rotational speed 
chosen. When the flow is assumed to follow the mean blade angle, an 
instantaneous change in flow angle and velocity occurs which results in 
an instantaneous rise in the absolute tangential velocity. 
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Unlike the flow in a stationary nozzle, the location of the points 
of minimum pressure, maximum velocity, and maximum Mach number in a 
rotating channel are not coincident and occur upstream of the geometric 
throat. The dashed lines in figure 3 connects the points of minimum and 

maximum values . 

The variation of the locations of the flow min i mum s and maximums 
is apparent from an inspection of the basic differential equations (5) 
to (7) from which the following equations axe obtained for zero slopes 
of p, V, and M in isentropic flow: 


Minimum pressure 



Maximum velocity 



Maximum Mach number 



(Ha) 


(lib) 


(11c) 


In the development of equations (ll), use was made of the equation 



which was derived from the general energy equation with the ratio of 
temperatures related to the ratio of the velocities of sound. 
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The locations of mini mum pressure and maximum velocity therefore 
apparently vary with the local flow Mach number as well as the channel 
geometry and rotor Mach number, hut the maxi mum local Mach number loca- 
tion is independent of flow. At the critical value of flow Mach num- 
ber (M = 1) , the locations of the pressure and velocity become coin- 
cident at the critical section. 

The location of the critical radius (M — 1) varies with both the 
rotor Mach number and the rate of convergence of the channel. With the 
use of a parabolic variation of area, several rates of convergence 
between the inlet and throat sections of the illustrative channel were 
obtained as the ratio of throat area to inlet area at constant throat 
area was varied. The effect of the rate of convergence and rotor Mach 
number on the location of the critical radius ratio is shown in figure 4. 
The curves illustrate the general consideration that because the con- 
traction and energy effects are in opposition, the largest displacements 
are to be expected for low rates of convergence and for high values of 
rotor Mach number. 

Maximum mass flow. - The maximum mass flow that a cross section of 
the channel is capable of passing under isentropic conditions is found 
by differentiating the expression for mass flow with respect to pres- 
sure ratio at the section and setting the derivative equal to zero. 

From this development (appendix B) it is found that the maximum mass 
flow occurs at the critical condition where the flow Mach number is 1 
and the pressure ratio is given by 



( 12 ) 


For radial flow with rotation, the critical pressure ratio increased 
with increasing impeller-tip Mach number. With no rotation, equa- 
tion (12) reduces to the well-known value for critical pressure ratio 
of 0.528 for energy-conservative systems (y = 1.40). The variation of 
critical pressure ratio with impeller -tip Mach number is shown in fig- 
ure 5 for the illustrative channel. 


The increase in critical density corresponding to the increase in 
critical pressure (fig. 5) indicates that the channel is capable of 
passing a greater mass of fluid as the rotor Mach number is increased. 
The variation of maximum mass flow as described by the flow function ^ 
is given by (appendix B) 
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The relation for the illustrative channel is shown in figure 6. The 
increase in maximum mass flow results not only from the increase in 
density but also from the increase in critical area as the critical sec- 
tion moves upstream from the throat. The curve of the ma.y-1 mim mass-flow 
function (fig. 6) is actually the locus of the mayimiTm points of the 
curves of variation of flow function with pressure ratio at various 
rotor Mach numbers. The general value of i|r as obtained from equa- 
tion (B7) of appendix B is 



The complete mass -flow spectrum for the illustrative channel up to an 
impeller tip Mach number of 1.8 is plotted in figure 7. 

Supersonic flow . - For values of outlet static pressure lower than 
that required for continuous isentropic compression at the critical 
mass flow, supersonic velocities are obtained in the passage down- 
stream of the critical section. As in the flow through stationary 
nozzles, the fluid expands isentropically to some point downstream of the 
critical section, passes through a normal compression shock, and is then 
further compressed isentropically to the prevailing outlet pressure. 

The location and strength of the shock are determined by the magnitude 
of the outlet pressure. The variation of static-pressure ratio and 
shock location in the supersonic -flow region of the illustrative chan- 
nel is shown in figure 8 for several supercritical outlet pressures. 

The dashed lines are curves of constant entropy corresponding to the 
various outlet pressures. The shock location for a given supercritical 
outlet pressure was determined by comparing the static -pressure ratio 
from the constant-entropy curve for that supercritical pressure to the 
isentropic expansion curve with the static-pressure ratio that would 
occur across a normal shock located at the particular radial position. 

An examination of the equations of motion and the entropy relation 
(as was done in reference 8) for flow with energy transfer due to 
radius change revealed that shocks are both mechanically and thermo- 
dynamically possible in the region between the critical section and the 
throat. Aside from the difference in location between the geometric and 
aerodynamic throats , all conventional theoretical supersonic flow and 
normal shock relations are valid for one -dimensional radial impeller 
flow. In the real impeller flow with finite shock thickness, the energy 
transfer may have some effect on shock stability, boundary- layer inter- 
action, and flow separation. 

Losses . - In the actual flow in channels of radial impellers, the 
existence of internal losses arising from friction, turbulence, and 
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mixing will modify to some extent the general picture of the flow 
behavior obtained in the previous sections. From equations (5) to (7), 
for adiabatic flow the effect of losses (by means of entropy and static- 
temperature increases) accelerated the flow along the channel, or in 
other words, losses can generally be considered to have the effect of a 
reduction in the flow area. As such, the critical radius, as well as 
the locations of minimum pressure and maximum velocity, will be dis- 
placed downstream nearer the geometric throat. 

The critical condition in the presence of losses continues to be 
characterized by the attainment of a relative Mach number of 1 as indi- 
cated by the singularity in equations (5) to (7) at M = 1. In the 
presence of losses, the pressure ratio at which sonic velocity occurs 
will be reduced and a corresponding reduction in max i mum mass flow 
will occur. The variation of maximum mass flow with critical pressure 
ratio (known as the choking line) for radial impeller flow expressed as 
the ratio of the mass-flow function at any irreversible condition to 
the ideal isentropic mass-flow function was found (appendix C) to be 



The maximum mass flow is, therefore, a function also of the location 
of the critical section and its variation with the losses. 


A relatively large effect of losses on the general flow behavior 
occurs when the channel flow involves a breakaway of the boundary layer 
or the formation of eddy regions. In such cases, the geometric-area 
variation of the channel can no longer be considered as the actual area 
governing the flow. For a given impeller channel, the variation of the 
actual flow area can be related to the principal flow variables from 
combination of the equation of continuity (equation (2)) and 

2 _ TP _ T dp 

" p “ ndp 

where n is the polytropic exponent, as 

dA dW 1 _ V_dV ( 16 ) 

A ’ " 4V a v 

P 


At this point, the concept of an energy number is introduced. The 
energy number is defined so that its square is equal to minus the ratio 
of the rate of change along the flow path of kinetic energy to the rate 
of change of static enthalpy, or 
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Equation (16) 


E 


2 


V dV 
dp 

P 


can therefore be expressed as 


dA dp / v E 2 \ 
A “ TP \n " ^2/ 


and, when the minimum or maximum of the flow area develops, the condi- 
tion will be characterized by 

M^^E 2 (17) 

Inasmuch as in most cases the effect on the average flow variation 
of the area reduction (which was due to breakaway) will be greater 
than the entropy increase accompanying the separation, the condition 
can be represented very closely by the simple relation 

M 2 - E 2 (18) 

Thus an intersection of the curves of the squares of Mach number and 
energy number along the channel should reveal the presence and location 
of a minimum or maximum of the actual flow area. 


EXPERIMENTAL ANALYSIS 

The use of one -dimensional theory in the investigation of the mean 
flow in channels of centrifugal impellers is illustrated by an analysis 
of the flow in the radial-inlet impeller of reference 5. 


Apparatus 

A photograph of the radial-inlet impeller is shown in figure 9. 

The direction of the flow at the entrance to and along the impeller was 
confined predominantly in a radial plane. The impeller contained 
18 blades and had an outer diameter of 14 inches. The impeller inves- 
tigation was conducted in a variable-component test setup with a 
34-inch vaneless diffuser at tip speeds from 800 to 1300 feet per sec- 
ond. In addition to standard performance instrumentation, static - 
pressure taps were installed along the impeller stationary front shroud. 
A complete description of the impeller design and experimental setup 
are given in reference 5. 
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Calculations 

The average state of the air at points along the channel of the 
impeller were calculated for the various tip speeds from the measure- 
ments of static pressure taken along the stationary-front shroud. 

Because of the essentially radial nature of the impeller configuration, 
the variation of pressure across the depth of the passage from front to 
rear shroud should he negligible compared with the variations along the 
streamline (as was indicated in reference 3 for a geometrically similar 
impeller) . The pressure readings along the stationary -front shroud 
therefore closely represent the average (one -dimensional) static pres- 
sure in the channel at each radial position. Other assumptions made in 
the calculations were: (l) The direction of the flow followed the cur- 

vature of the blades at all points, (2) the impeller channel was flowing 
full at all times, and (3) the flow was isentropic. The detailed pro- 
cedure involved in the calculating of the principal flow variables from 
the continuity condition and energy equation (similar to equation (10)) 
are given in reference 5. 


Analysis of Flow 

General variations . - Illustrative examples of the measured varia- 
tion of static pressure along the impeller radius and the calculated 
variations of velocity and Mach number are shown in figure 10 for four 
weight flows at an impeller -tip Mach number of 0.878. Although the 
measured static-pressure variation with radius was not precisely defined 
in the inlet region of the impeller because of the few tap locations 
used in the investigation presented in reference 5, subsequent unpub- 
lished data from the same impeller with a considerably larger number of 
pressure taps along the front shroud have substantiated the faired var- 
iations used for the original data. 

The plots of static pressure, velocity, and flow Mach number 
(fig. 10) reveal a general similarity to the flow in the illustrative 
convergent -divergent channel in figure 3. The effect of the formation 
of a throat in the inlet region due to increasing relative approach 
area and leading edge separation at high weight flows is apparent. In 
the maximum flow range, the pressure drop became quite pronounced, 
resulting in local Mach numbers becoming and exceeding 1. After the 
maximum weight flow was attained, further reductions in discharge pres- 
sure increase the magnitude of the supersonic velocities and accompanying 
shocks and shift the shock regions further downstream from the throat. 

As a result of these effects, a considerable decrease in efficiency was 
obtained. 
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Effective flow area . - In a previous section a minimum or maximum 
of the flow area is shown to be characterized by the equality of the 
squares of the local Mach number M 2 and the energy number E 2 . For 
the geometric -area variation of the test impeller (fig. ll) , if the 
effective channel in the impeller -inlet region is to be considered non- 
rigid and varying with weight flow such that a throat forms at the 
higher weight flows, then the qualitative radial variation of the prin- 
cipal flow variables over the complete weight-flow range at constant 
tip speed will appear as illustrated in figure 12. The corresponding 
qualitative variations of Mach number and energy number will be as shown 
in figure 13. One intersection should be anticipated in the outlet 
region where the maximum geometric area occurs, and in the high flow 
range, a second intersection should be expected in the inlet region as 
the throat develops. 

Values of Mach number and energy number were calculated for the 
test impeller based on the ideal geometric area of the impeller channel 
and are shown in figure 14 for the test range of flow at an impeller- 
tip Mach number of 0.878. The calculated values of M 2 and E 2 in 
general follow the qualitative variations of figure 13 rather closely. 

The entrance -region intersection, as expected, indicates the existence 
of a minimum of the flow area. The downstream shift of the minimum of 
the flow area with an increase in flow is also apparent. Similar var- 
iations of E 2 and M 2 were obtained for the runs of an impeller-tip 
Mach number of 0.965. 

Calculations of the effective flow area which was required to pass 
the measured maximum weight flow and to produce the measured pressure- 
ratio variation were made for the accelerating flow region at the inlet 
for several impeller -tip Mach numbers. The effective flow area was 
determined from equation (10) in which the area A was introduced by 
means of the continuity substitution for the relative velocity V. A 
polytropic efficiency of 0.95 was assumed for the pressure-density 
relation. The calculated flow area in the accelerating flow region at 
maximum weight flow at an impeller-tip Mach number of 0.878 is shown in 
figure 11 in comparison with the geometric area and relative approach 
area in order to illustrate the nature of the effective area contraction. 

Maximum weight flow. - At a given tip speed, the magnitude of the 
maximum weight. flow resulting from the attainment of the critical con- 
ditions within the impeller passage depends on the location of the 
critical section and efficiency of the flow. As a first approach, ideal 
critical conditions were assumed. The location of the critical radius 
was obtained for each tip-speed run by comparing the calculated ideal- 
critical pressure ratio as given by equation (12) with the measured 
static -pressure ratios (fig* 15 for impeller-tip Mach numbers of 0*878 
and 1*058) * The maximum corrected weight flow at each tip speed was 
obtained from equation (13) and the calculated effective flow area 
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(as in fig. 11) at the critical radius. Maximum weight flow plotted 
against impeller-tip Mach number is shown in figure 16 in comparison 
with the measured values of reference 5. As expected, the calculated 
maximum weight flow agreed very closely with the experimental values. 

Values of maximum weight flow predicted by the conventional prac- 
tice of assuming the choking to occur at the impeller inlet (vane lead- 
ing edge) are also shown in figure 16. The assumption of inlet choking 
apparently is inaccurate for predicting maximum weight flow for this 
type of impeller, and an improved method based on choking at the effec- 
tive throat within the passage is desirable. 

In plotting the calculated effective flow area, the effective flow 
area at the critical radius was observed to be very nearly equal to the 
geometric area at all tip speeds (for example fig. ll) . This was 
believed to be due to the consideration that, as the flow area varied 
from the relative approach area upstream of the inlet (greater than the 
geometric area) to the throat area in the separated region downstream 
of the inlet (smaller than the geometric area), there is some point 
upstream of the throat where the effective flow area is equal to the 
geometric area. With the critical section occurring upstream of the 
throat, the displacement from the throat was of sufficient magnitude to 
result in the near equality of the areas. Thus, for a given impeller, 
if the location of the critical radius can be approximated (from an 
inspection of the geometric -area variation and the relative approach 
areas) , the ideal weight flow based on geometric area at the critical 
section may provide a method of predicting maximum weight flow that is 
more accurate than the inlet-choking estimate. For the radial-inlet 
impeller of reference 5, ideal maximum weight flow based on geometric 
area at the critical section is shown by the dashed line in figure 16. 

Losses and efficiency . - The measured variation of over -all effi- 
ciency for the radial-inlet impeller installation of reference 5 as 
plotted against relative approach area and angle of attack is shown 
over a range of impeller -tip Mach numbers in figure 17. Peak efficiency 
for all impeller -tip Mach numbers occurred at values of relative 
approach area from about zero to 10 percent greater than the geometric 
area of the channel at the blade inlet (including blade thickness) 
corresponding to positive angles of attack from 3 to 6 degrees. Angle 
of attack is defined as the difference between the angle of the inlet 
air and the tangent to the blade mean line at the leading edge. The 
direction of the approaching fluid was calculated on the basis of inlet- 
depression tank stagnation conditions and the full area of the flow 
passage immediately upstream of the leading edge of the impeller blades. 
A 3 -percent velocity- increase correction for shroud boundary- layer 
blockage and total-pressure loss was included. 
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The over -all efficiency variation shown in figure 17 is a result of 
the combined effects of the various losses occurring in the flow. In 
general, heavy contributors to impeller losses can be identified as 
(l) friction losses arising from the action of the mean flow on the 
wetted surface of the impeller and vaneless diffuser passages, (2) 
losses due to the somewhat sudden enlargement or contraction of the rel- 
ative flow area across the blade inlet (referred to as inlet losses), 

(3) boundary-layer and secondary-flow losses due to the effects of 
blade-surface velocity gradients arising from the blade loading 
(referred to as blade -loading losses), and (4) shock losses due to 
supersonic velocities in the choked-flow range of operation. Friction 
losses depend principally on the square of the mean flow velocity and 
will, therefore, increase rapidly with increasing weight flow. Blade- 
loading losses, on the other hand, increase with decreasing weight flow 
as there is less kinetic energy available in the mean flow to overcome 
the adverse pressure gradients and boundary-layer growth imposed by the 
blade circulation. (For a radial discharge impeller with little slip, 
the circulation is approximately constant at a fixed tip speed.) 

The relative magnitude of the inlet losses depends largely on 
whether the relative inlet area is contracting or expanding across the 
blade inlet. For low weight flows the flow area is enlarged (fig. 1), 
and from fluid mechanics for pipe flow (reference 9 or 10), the loss 
for sudden enlargement is given by (neglecting density changes) 


In this loss analysis, subscript 1 refers to conditions ahead of the 
blade leading edge and subscript 2 refers to conditions in the channel 
imme diately downstreami of the leading edge. For contracting areas at 
high weight flows, 



where the factor K is a function primarily of area ratio. For values 
of A 2 /Ap to about 0.5, the loss can be given approximately by 


Thus, as flow is decreased from the point of equal inlet area (minimum 
entrance loss) , although the area ratio A]_/a 2 decreases in equa- 
tion (19), the relative inlet velocity V]_ also decreases and the 



(19) 



( 20 ) 
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entrance losses will remain very small. In contrast, however, as flow 
is increased past the equal-area point, Ag/A]_ decreases and 

increases (V 2 approaches the critical value) and a sharp rise in inlet 
loss occurs. 

At a given impeller -tip speed, the composite loss picture in fig- 
ure 17 is as follows: On the low weight-flow side of the equal-area 

point, friction losses decrease gradually, entrance losses remain very 
small, blade loading losses tend to increase and a gradual decline in 
efficiency to the surge point occurs. On the high weight-flow side, 
blade loading losses gradually decrease, but friction and inlet losses 
rise sharply to produce a rapid rate of decrease in efficiency to the 
choked flow point. The large declines in efficiency at the choked flows 
at the higher tip speeds are due to the shock losses of the super- 
critical passage flow. 

The general decrease in efficiency as impeller-tip Mach number is 
increased (fig. 17) is a reflection of the increased friction and 
blade -loading losses resulting from the higher values of circulation 
and mean velocity at these wheel speeds. The decrease in range of 
operation at the higher tip speeds can also be attributed to the higher 
velocities and loadings. In the high flow region of operation, as the 
relative inlet Mach number is increased because of the higher wheel 
speeds, the area reduction necessary to cause a critical contraction 
ratio becomes progressively smaller. 

At the choked flow point (local sonic velocity in the blade inlet 
region) the over -all efficiency is substantially greater at the higher 
wheel speeds than it is at the lower wheel speeds. The reduced over- 
all efficiency at critical flow at low wheel speeds cannot be attri- 
buted to increased diffuser, friction, or blade-loading losses, because 
(l) impeller channel velocities are not substantially different for the 
various wheel speeds at the critical condition", (2) impeller blade 
loadings remain light at the lower wheel speeds; and (3) vaneless dif- 
fuser efficiency gives evidence of little change with wheel speed (ref- 
erences 11 and 12) . A principal difference between the high- and low- 
speed choking points however, is the flow-area change across the blade 
inlet. According to equation (20), with the velocities downstream of 
the area contraction approximately the same (local critical velocity), 
a greater inlet loss occurs at the lower wheel speeds where the inlet 
area change is the greatest. A similar condition of large inlet losses 
in the high flow range of a geometrically similar impeller at low wheel 
speeds was found in references 3 and 4. 

In view of the indicated large inlet losses at negative angles of 
attack, improved performance and range of operation may be obtained if 
blade leading-edge shapes are designed for a more efficient accomoda- 
tion of the inlet -flow variation. 
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SUMMARY OF RESULTS 

From a theoretical analysis of the ideal one -dimensional flow in 
9-11 illustrative rotating radial impeller channel with a convergent- 
divergent area, it was found that: 

1. The behavior of the flow (variations of pressure, velocity, Mach 
number, etc.) along the radius of the rotating channel was generally 
similar in trend to that in a geometrically similar stationary nozzle. 

2. Critical (sonic) conditions occurred upstream of the geometric 
throat, with the displacement from the throat increasing with increasing 
impeller wheel speed and with decreasing rate of area contraction. 

3. The effect of losses on the variation of flow along the channel 
was similar to the effect of a reduction of the flow area. 

A further one -dimensional analysis of the flow in an experimental 
radial-inlet impeller containing static-pressure taps along the sta- 
tionary front shroud indicated that; 

1. The behavior of the mean flow along the impeller passage was 
generally similar to that of flow along a rotating radial channel whose 
effective flow area at the inlet varies with the operating point. 

2. A critical condition occurred in the inlet region of the impeller 
at high negative angles of attack where a throat is formed because of 
flow separation and readjustment downstream of the leading edge. 

3. Calculated maximum weight flows based on conditions at the 
critical radius compared favorably with experimental maximum weight 
flows over the wheel speed range of the impeller. 

4. Impeller -inlet losses were indicated to be large at the higher 
weight flows. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
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APPENDIX A 


LOGARITHMIC DIFFERENTIAL EQUATIONS 

The variations of pressure, relative velocity, and relative Mach 
n umb er along the radius of a channel can be obtained from equations (1) 
to (4) in the logarithmic differential form as follows. 

Velocity. - From equation (4) , 


dt = i d / £ 
R l P 


Substituting in equation (l) yields 


VdV ♦ £ d (p) - d « ♦ 


Differentiating p/p, with R = Cp - c v , gives 


VdV + 


r-i 


4p _ / P\ *£ 

P I P7 P 


= dQ + d 


ao 2 r 2 


From continuity, equation (2), 


(Al) 


(A2) 


dp _ dA dV 

p A ' V 


(A3) 


Substituting for and ~~ in equation (A2) from equations (A3) and 

(3) yields, with a 2 = 


1 V 2 dV a 2 dV 1 d /a> 2 r 2 \ . X « * _dA 

7^1 — + 7=uT =dQ_ Fi d I — / + 7=1 ^ A 


r-1 9 . 

Multiplying through by — o— , using M = — o and condensing give 
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dV 

V 


1 

l-l£ 



u> 2 r 2 


-ydF 




Press-ore, - From eq-oation (3), using 
yields 


2 

a 


IP 

P 


and equation 


- 1 /JL frl ( (X)2r2 \ 

p " i-rf U 2 L \ 2 / 


[(r-l)J^ + 1] dF - (r-l)M 2 dQ 


+ 


Mach number. - From the definition of Mach number M = V/a 
he shown that 


dM dV 1 da^ 

M V " 2 q 2 

ct 


From the definition of sonic velocity. 



Thus, equation (A4) can he expressed as 


dM dV 1 /dp dp\ 
M V - 2 yp ” p / 


Substituting for 


dp 

P 


from equation (A3) gives 


dM _ 1 /dV _ dp dA\ 

M 2\V~p~A/ 

Substituting for —■ and — in equation (A7) from equation (6) 
then yields after combining terms 


( 6 ) 

( 6 ) 



it can 


(A4) 

(A5) 

(A6) 

(A7) 
and (5) 
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cLM 

M 


1 

2(l-Mp) 



(1+r) CL 



r [2 + ( r -i) wP] dF - (r-i)(i + tmP) <iQ 


+ [2 + (r-i) x) 


( 7 ) 
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APPENDIX B 


IDEAL MAXIMUM MASS FLOW 


From continuity, the ideal mass flow at any cross section of a 
channel for one- dimensional flow is given by 


W = pVA 


(Bl) 


From the energy equation 


V = 


2c p JT 0 + <i> 2 r 2 - 


From the isentropic relation for ideal flow 

1 

P = Pn I 


o p 


(B2) 


(B3) 


Thus, 


r p 
r-i p 


r o /p 

T - 1 Pq V P o 


hi 

r 


(B4) 


Now, inasmuch as 


o 2 


equation (B4) can he expressed as 


2 

_r_P T 

r-i p r-i Ip 0 


(B5) 


Furthermore, from the perfect gas law, 


P a 2 
c JT = JT _o o 

P 0 r-1 p r-1 

K o 


(B6) 
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Substituting equations (B2), (B3) , (B5) , and (B6) into equation (Bl) 
results in 


W = Po 



A 


2 

T- 1 


a o 


a> 2 r n 



2a 0 2 

r-i 


r-ii 1 



from which is obtained 


W = 



a 0 P 0 




(B7) 


Diff erentiating W with respect to pressure ratio at a given cross- 
sectional area yields 



Inasmuch as the denominator is a nonzero quantity (see equation (B7)), 
it is sufficient to set only the term between braces in the numerator 
equation zero in order to maximize equation (B7 ) , or 


2-r 1 

2 

r 

from which the pressure ratio at maximum flow is obtained as 




2 . r-i 

j+i r+i 



(B12) 


The velocity at the maximum flow condition is found by expressing 
equation (B2) in terms of the local velocity of sound 
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V m =/ 


1 2 

r-l / r \ 2 M 2 a m 2 

Jr - 1 a ° 

L 1 + - y m % J 


Then, from equation (B8) 


* m 


a ° 


y+1 /p 
2 VP,. 


r-i 

r a 


m 


Substituting the adiabatic relation 

r-i 


if 

in a o 


* 0 / 


m 


(B9) 


into equation (B9) and rearranging give 


V - (l+± o 2 _ -L. * 

m \r-i ^ r-i ^ 


from which is obtained the critical condition 


^m - a c 

Thus the maximum mass flow occurs at the critical condition and all m 
subscripts can be replaced by c subscripts. 

The maximum mass flow is given by 


W — Pc a c 


which becomes, after the introduction of the isentropic and sonic 
relations 


Wm - po 3-o Aq 



(BIO) 
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Substituting for the critical pressure ratio from equation (12) 


W 


m 


p o a o^c 


= t m = 





y+1 

ZCFD 


(13) 
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APPENDIX C 


MASS -FLOW CHOKING LINE FOR FLOW WITH LOSSES 

Maximum mass flow in a rotating channel in the presence of losses 
can be obtained as follows. 


From continuity, at the critical condition 

t IT! 

W c ~ Pc a c ^c 

With the use of the sonic relation and equation of state. 


(Cl) 


K = 


P’V?K 


A* 


(C2) 


From the energy equation for relative flow, 

a c 2 o> 2 r 2 

c P Jt c + — = V^o + 


or 


Cp Jt C + 


rst, 


c_ n , ^ 2 r c 2 

2 “ + 2 


from which 


or 


2 / v- i (X ^ r 2 

tc = f+1 T ° ( 1 + ~ “ 2" 


s 7+1 T o 


1 + 


r-i/£c 
r T 




(C3) 


Substituting equation (C3) into equation (C2) gives 


P 1 
*0 


a/tRTq A c P q 


W c 3 


RT, 


r+i 


1 + r iK|)V 
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and 


W 

c 


Po -^c 


= -Or’ r= 
Y c 



(C4) 


The isentropic critical mass-flow function is obtained from equa- 
tion (BIO) as 


r-i 



and from equation (12) for critical pressure ratio 


(C5) 



(C6) 


The ratio of actual to isentropic critical mass-flow function or choking 
line thus becomes 



(15) 
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Rotation 



(a) Approach area less than blade area 
(low weight flow). 



(b) Approach area equal to blade area 
(medium weight flow) . 



(c) Approach area greater than blade 
area (high weight flow) . 


Figure 1. - Variation of relative approach area over 
weight-flow range of impeller operation. 
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(a) Flow area variation. 




(to) Equivalent circular nozzle. 

Figure 2. - Illustrative impeller channel for theoretical 

flow analysis. 
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Figure 2. - Concluded. 


Illustrative impeller channel for theoretical flow analysis. 


Static -pressure ratio, P/P 
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(a) Static -pressure ratio. 


Figure 3. - Variation of flow along illustrative impeller channel from zero 
to critical flow. Impeller-tip Mach number My, 1.0. 


Velocity, V, ft /sec 
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(t) Velocity. 

Figure 3. - Continued. Variation of flow along illustrative impeller channel 
from zero to critical flow. Impeller-tip Mach number Mr, 1.0. 


Mach number 


oo 

-F=" 



( c ) Mach number . 

Figure 3. - Concluded. Variation of flow along illustrative impeller channel 
from zero to critical flow. Impeller-tip Mach number M^ 1.0. 
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Critical radius ratio, (r/r T ) 



Figure 4. - Variation of critical radius ratio with impeller-tip Mach number for various throat-area ratios ^ 

(constant throat area). vji 
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0 .4 .8 1.2 1.6 2.0 

Impeller-tip Mach number, Mij 

Figure 5. - Variation of critical pressure ratio with impeller- 
tip Mach number for illustrative impeller channel. 



. 1— — _ — j 1 

0 .4 .8 1.2 1.6 2.0 

Impeller-tip Mach number, Mrp 


Figure 6. - Variation of maximum mass-flow function with impeller- 
tip Mach number for illustrative impeller channel. 


Mas s -flow function 
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Figure 7. - Variation of mass-flow function with pressure ratio at critical 
radius for illustrative impeller channel. 


Static -pres sure ratio, p/P ( 
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Radius ratio, r/r^, 


Figure 8. - Variation of static-pressure ratio along illustrative impeller 
channel for several supercritical out let -pres sure ratios. Impeller-tip 
Mach number Mfj>, 1.0. 
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Figure 9. - Experimental radial- inlet impeller (reference 5). 


Mach number, M Velocity, V, ft/eec Static -pressure ratio, p/P ( 
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(a) Static -pressure ratio. 



Radius, r, in. 
(c) Mach number . 


Figure 10. - Measured variations of shroud static 
pressure and calculated variations of velocity 
and Mach number along channel of experimental 
radial- inlet impeller. Impeller-tip Mach num- 
ber, 0.878. 
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12 i ^ ^ ^ 1 ^ ^ 1 1 

3 4 5 6 7 


Radius, r, in. 

Figure 11. - Variation of flow area along radius of 
experimental radial-inlet impeller at maximum 
weight flow. Impeller-tip Mach number Mr, 0.878. 
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Figure 12. - Qualitative variation of one-dimensional effective flow area A, static pressure p, velocity V, and potential energy 
along channel of experimental radial-inlet impeller as flow Is increased from minimum to maximum at constant impeller-tip speed. 
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Figure 13. - Qualitative one-dimensional variation of square of Mach number and energy number E* along chan] 

radial-inlet impeller as flow is increased from minimum to maximum at constant impeller-tip speed 
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Square of Mach or energy number 
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(a) Weight flow, 
2.774 pounds 
per second. 



4.370 pounds 
per second. 



2 4 6 8 

r, in. 

(c) Weight flow, 


4.795 pounds 
per second. 



(d) Weight flow, 
5.309 pounds 
per second. 


Figure 14. - Calculated variations of square of Mach number M 2 and energy number E 2 along channel of experimental radial-inlet 
impeller as flow is Increased from minimum to maximum. Impeller-tip Mach number Mip, 0.878. 
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Square of Mach number or energy number 


2.0 



(e) 


Weight flow, 


5.444 pounds 


per second. 



2 4 6 8 

Eadius, r, in. 

(f) Weight flow, 5.501 pounds 


per second maximum flow. 



2 4 6 8 


(g) Weight flow, 5.428 pounds 

per second, supercritical P^/Pq» 


Figure 14. - Concluded. Calculated variations of square of Mach number M 2 and energy number E 2 along channel of 
experimental radial-inlet impeller as flow is increased from minimum to maximum. Impeller -tip Mach number 
Mr, 0.878. 
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Figure 15. - Variation of static-pressure ratio along radius of experimental radial-inlet impeller. 
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Adiabatic efficiency, rj a ^ Maximum weight flow, lb /sec 
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Figure 16. - Variation of maximum weight flow with impeller-tip Mach number 
for experimental radial-inlet impeller. 



Relative approach area at inlet A-p sq in. 


16 12 8 4 0 -4 -8 -12 

Angle of attack^ deg 

Figure 17. - Variation of over-all adiabatic efficiency with relative 
approach area and angle of attack at inlet of experimental radial- 
inlet impeller. 
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